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Abstract 

We investigate the constraints on the flavour violating parameters from the 
decay B — ► X s ~f, taking into account the interplay of the various sources of 
flavour violation in the unconstrained MSSM. We present a systematic leading 
logarithmic QCD analysis of these model-independent constraints, including 
contributions from gluinos, neutralinos, charginos, charged Higgs bosons and 
interferences between them. We show that two simple combinations of ele- 
ments of the down squark mass matrix are stringently bounded over large 
parts of the parameter space where only weak assumptions on the hierarchi- 
cal structure of the squark mass matrices are made. We also briefly analyse 
up to which values SUSY contributions, compatible with B — > X s j, can en- 
hance the Wilson coefficient Cs(mw), which plays an important role in the 
phenomenology of charmless hadronic B decays. 
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I. INTRODUCTION 



Today supersymmetric models are given priority in the search for new physics beyond 
the standard model (SM). This is primarily suggested by theoretical arguments related to 
the hierarchy problem. Supersymmetry eliminates the sensitivity for the highest scale in the 
theory and, thus, stabilizes the low energy theory. 

Flavour changing neutral current (FCNC) processes provide crucial guidelines for super- 
symmetry model building. In the so-called unconstrained minimal supersymmetric stan- 
dard model (uMSSM) there are new sources for FCNC transitions. Besides the Cabibbo- 
Kobayashi-Maskawa (CKM)-induced contributions, there are generic supersymmetric con- 
tributions induced by flavour mixing in the squark mass matrices. The structure of the 
uMSSM does not explain the suppression of FCNC processes, which is observed in experi- 
ments; this is the crucial point of the well-known supersymmetric flavour problem. Within 
the framework of the MSSM there are at present three favoured supersymmetric models that 
solve the supersymmetric flavour problem by a specific mechanism through which the sector 
of supersymmetry breaking communicates with the sector accessible to experiments: in the 
minimal supergravity model (mSUGRA) |]J, supergravity is the corresponding mediator; 
in the other two models, this is achieved by gauge interactions and by anomalies [[3]. 
Furthermore, there are other classes of models in which the flavour problem is solved by 
particular flavour symmetries Q]. 

Flavour violation thus originates from the interplay between the dynamics of flavour 
and the mechanism of supersymmetry breaking. FCNC processes therefore yield impor- 
tant (indirect) information on the construction of supersymmetric extensions of the SM 
and can contribute to the question of which mechanism ultimately breaks supersymmetry. 
The experimental measurements of the rates for these processes, or the upper limits set 
on them, impose in general a reduction of the large number and size of parameters in the 
soft supersymmetry-breaking terms present in these models. Among these processes, those 
involving transitions between first- and second-generation quarks, namely FCNC processes 
in the K system, are considered as the most formidable tools to shape viable supersym- 
metric flavour models. Moreover, the tight experimental bounds on some flavour-diagonal 
transitions, such as the electric dipole moment of the electron and of the neutron, as well as 
g — 2, help constraining soft terms inducing chirality violations. 

Among neutral flavour transitions involving the third generation, the rare decay B — > 
X s 7 is at present the most important one, as it is the only inclusive mode which is already 
measured pj. The theoretical SM prediction, up to next-to-leading logarithmic (NLL) pre- 
cision || for its branching ratio, is in agreement with the experimental data. Although the 
experimental errors are still rather large, this decay mode already allows for theoretically 
clean and rather stringent constraints on the parameter space of various extensions of the 
SM (see for example [0]). 

Once more precise data from the B factories are available, this decay will undoubtedly 
gain efficiency in selecting the viable regions of the parameter space in the above classes of 
models; it may also help discriminating among the models that will be proposed by then. 
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In this paper we present a model-independent analysis of the decay B — > X s 7, based on a 
LL-QCD calculation, where contributions from IT^-bosons, charged Higgs bosons, charginos, 
neutralinos and gluinos are systematically included. Former analyses in the unconstrained 
MSSM neglected QCD corrections and only used the gluino contribution to saturate the 
experimental bounds. Technically, the so-called mass insertion approximation (MIA) was 
used where the off-diagonal elements of the squark mass matrices are taken to be small and 
their higher powers neglected. As a consequence of this single insertion approximation, no 
correlations between different sources of flavour violation were taken into account. In this 
way, one arrived at 'order-of-magnitude bounds' on the soft parameters PHTT|. In []Efl , the 
sensitivity of the bounds on the down squark mass matrix to radiative QCD corrections was 
analysed including the SM and the gluino contributions. The aim of the present paper is to 
extend this analysis to include the contributions from charged Higgs bosons, charginos and 
neutralinos and their interference effects, and even more important, the effects that result 
when several flavour violating parameters, i.e. several off-diagonal elements in the squark 
mass matrices, are switched on simultaneously. We anticipate that two simple combinations 
of matrix elements of the down squark mass matrix remain rather stringently bounded 
over large parts of the parameter space, in a general scenario where only relatively weak 
assumptions on pattern of the squark mass matrices are made. 

Since there are different contributions to this decay, with different numerical impact on 
its rate, some of these flavour- violating terms may turn out to be poorly constrained. Thus, 
given the generality of such a calculation, it is convenient to rely on the mass eigenstate 
formalism, which remains valid even when some of the intergenerational mixing elements 
are large, and not to use the approximate mass insertion method, where the off-diagonal 
squark mass matrix elements are taken to be small and their higher powers neglected. In 
the latter approach the reliability of the approximation can only be checked a posteriori. 

Finally, we note that the off-diagonal elements of the squark mass matrices can get 
constraints on completely different grounds, namely from the requirement of the absence 
of charge and colour breaking minima as well as from directions in the scalar potential 



which are unbounded from below (see ||T3| for a more detailed discussion). However, these 
bounds correspond to sufficient, but not necessary conditions for the stability of the standard 
vacuum, because it is possible that we live in a metastable vacuum, whose lifetime is longer 
than the age of the universe [fPl ]. 



The paper is organized as follows: In section 2, we discuss the framework for the calcula- 
tion of the branching ratio for B — > X s / y. In section 3, we briefly recall the sources of flavour 
violation, encoded in the squark mass matrices. In section 4, we present the phenomeno- 
logical analysis on the bounds on the flavour violating parameters. In section 5, we briefly 
explore up to which values SUSY contributions, allowed by B — > X s 7, can enhance the Wil- 
son coefficient Cg(mjy), which plays an important role in the phenomenology of charmless 
hadronic B decays. In section 6 we give a short summary. In appendix Al, we state our 
conventions, while in appendix A2 we list the Wilson coefficients at the matching scale. 



II. FRAMEWORK FOR CALCULATING B -> X S J 
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A. Hamiltonians 



In the SM, rare 5-meson decays are induced by one-loop diagrams in which W bosons 
and up-type quarks propagate. The most important corrections to the decay amplitude 
for b — > S7 are due to exchanges of gluons, which give rise to powers of the factor L = 
\og(ml / rri^y) . It turns out that each of these logarithms is accompanied by at least one 
factor of a s . Since the two scales m b and m w are far apart, L is a large number and 
these terms need to be resummed: at the leading logarithmic (LL) order, powers of a s L 
are resummed; at the next-to- leading (NLL) order, also the terms of the form a s (a s L) N 
are systematically resummed. Thus, the contributions to the decay amplitude are classified 
according to 

(LL) : G F (a s L) N , (NLL): G F a s (a s L) N , (N — 0, 1, ...) (1) 

where Gf is the Fermi constant. The resummation of these corrections is usually achieved by 
making use of the formalism of effective Hamiltonians, combined with renormalization group 
techniques. The effective Hamiltonian Ti^f, obtained by integrating out the top-quark and 
the W boson, can be written as 

AG 

nZf = ~^ K A E CiMOiM . (2) 

The Wilson coefficients C{ contain all dependence on the heavy degrees of freedom, whereas 
the operators Oi depend on light fields only. The operators relevant to b — > 57 read 

1 = (s lp T a P L c) (&f»T a P L b) , 

02 = (s^Plc) (cYP L b) , 

3 = (s llx P L b)E q (qrq), 

4 = (s^T a P L b)E q (qrT a q), 

05 = (s llllvlp P L b) E q (qrYYq) , ( 3 ) 

O e = (s^ ulp T a P L b)E q (qrYYT a q), 
e 

16V 



°7 = 7^ m b(v)(sa^P R b)F, u , 



Os =-^m b (^)(sa^T a P R b)G; i/ . 

The matrices T a (a = 1, ...,8) are SU(3) colour generators and P LjR are left- and right- 
handed projection operators; e and g s denote the electromagnetic and the strong coupling 
constants, respectively. Note that the 6-quark mass is the relevant parameter that governs 
the chirality flip in the SM dipole operators O7 and 0$. All eight operators are of dimension 
six. We anticipate that this is in contrast with the dipole operators induced by gluinos, 
where the helicity flip can be generated by the gluino mass instead of the 6-quark mass, as 
we will see in more detail later. As a consequence, these dipole operators are effectively of 
dimension five. 

A consistent SM calculation for B — > X s j at LL (or NLL) precision requires three steps: 
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1) a matching calculation of the full SM theory with the effective theory at the scale 
H = Hw to order a° s (or a\) for the Wilson coefficients, where nw denotes a scale of 
order m w or m t ; 

2) a, renormalization group treatment of the Wilson coefficients using the anomalous- 
dimension matrix to order a] (or a 2 s )\ 

3) a calculation of the operator matrix elements at the scale // = /i& to order ot° s (or a]), 
where //& denotes a scale of order m&. 

In supersymmetric models there are further contributions to the FCNC processes studied 
in this paper, i.e. the exchange of charged Higgs bosons and up-type quarks; of charginos 
and up-type squarks; of neutralinos and down-type squarks; and of gluinos and down-type 
squark. They all lead to |A(5)| = |A(,S)| = 1 effective magnetic and chromomagnetic 
operators (of Cy-type, CVtype) and also to new four-quark operators. 

Taking into account operators up to dimension six only, the effects of charginos, neu- 
tralinos and charged Higgs bosons can be matched onto the usual SM magnetic and chro- 
momagnetic operators and Og and their primed counterparts 



°' 7 = lib mbifl) ^ Pl ^ F ^ ' = 1^2 m M (s^T a P L b) G% . (4) 



We would like to stress that we do not work in the mSUGRA scenario. Therefore, a lot of 
the relations deduced in || do not hold anymore. The most important thing to notice is 
that in a general SUSY framework, there is no connection between the Yukawa couplings 
in the superpotential and the corresponding trilinear term in the soft potential. However, 
the contributions from charginos and charged Higgs bosons to the Wilson coefficients of the 
primed operators vanish also in the more general unconstrained MSSM if m s is put to zero 
(see eqs. (|A18 ) and (|A19 )). This implies that for physical values of m s the chargino- and the 



charged Higgs boson contributions to the primed operators are small in scenarios in which 
tan/3 does not take extreme values. The neutralino contributions to both, the primed and 
unprimed operators are also unimportant, because their Wilson coefficients involve those 
entries of the squark mixing matrices Tdl and Tor, which also induce gluino contributions; 
the latter, which are proportional to g 2 therefore dominate the neutralino contributions 
which are proportional to g\. We also found numerically that the neutralino contibutions 
are indeed rather inessential. 

The fact that the operators generated by charged Higgs bosons, charginos and neutralinos 
can be matched onto the SM operators and their primed counterparts implies that the terms 
that get resummed at LL, show the same pattern as those listed in eq. ([]]); the Fermi 
constant appearing there is obviously replaced by a specific supersymmetric parameter for 
the chargino and neutralino contributions. 

Matters are somewhat different for the gluino contribution Ti. 9 e ff^ as worked out in detail 



in 



12[ . At the amplitude level, terms of the form 

(LL): a s (a s L) N } (NLL): a s a s (a s L) N , (iV = 0,l,...) 
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are resummed, respectively at the leading and next-to-leading order. While Ti-eff i s unam- 
biguous, it is a matter of convention whether the a s factors should be put into the definition 



of operators or into the Wilson coefficients. We follow the framework developed in ref. [12 



where the distribution of the a s factors was done in such a way that the anomalous dimen 



sion matrix systematically starts at order a s . We write the effective Hamiltonian Ti,^, ff in 



the form 

Kff = E c iMOiM + EE C !M°IM ■ ( 5 ) 

i i q 

The index q runs over all light quarks q = u,d, c, s, b. Among the operators contributing to 
the first part, there are dipole operators in which the chirality flip is induced by the 6-quark 
mass: 

lb . g = e g 2 {fi) m b {fi) {sa^P R b) , 0' 7bg = e g 2 {fi) m 6 (» (sa^P L b) , 

(6) 

%b - g = g 3 s ( f x)m b (fi)(sa^T a P R b)G« v , Q' m = <? s 3 (/i) m 6 (//) (sa^T a P L b) G% . 



As discussed in |12[ , there are also gluino- induced operators where the chirality violation is 
signalled by the charm quark mass (obtained by replacing m b (n) by m c (/i)) and operators 
where the chirality flip is induced by the gluino mass. The latter read 

7g , g = eg 2 M (sa^P R b) F, v , 0' 7M = eg 2 ^) (sa^P L b) F, v , 

7 

8g , g = g 3 s (fi) (sa^T a P R b) , & m = g 3 (fi) (sa^T a P L b) G% . 

At the LL- level, these operators could be absorbed into the operators given in eq. (d), when 
neglecting the small mixings effects from the gluino- induced four- Fermi operators with scalar 
or tensor Lorentz structure. However, as it is useful to separate the contributions where the 
chirality flip is induced by m g , we do not perform this absorption. Notice that the operators 
in eq. (0) have dimension five, while the other operators are of dimension six. We also stress 
that unlike the other supersymmetric contributions, the primed gluino-induced operators 
are not suppressed compared with the unprimed ones. This is in strong contrast with the 
mSUGRA scenario, where the primed operators are stronlgy suppressed ||. 

The contributions to the second part in eq. (|j) are given by four-quark operators with 
vector, scalar and tensor Lorentz structure. As shown in ref. 0], the scalar and tensor 



operators mix at one loop into the six- dimensional magnetic and chromomagnetic ones. 
Therefore, they have to be included in principle in a LL calculation. As mentioned above, 
these mixings are numerically small and therefore not very important in practice. 

For completeness we recall all Wilson coefficients at the matching scale fiw i n ap- 
pendix [A 2j . The anomalous dimension matrix of the SM operators Oi-Os and the evolution 
of the corresponding Wilson coefficients to the decay scale fi b are well known and can be 
found in ||. The evolution of the gluino-induced Wilson coefficients C^g is given in ref. [I2~|. 
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B. Branching ratio 



The decay width for B —> X s ^y to LL precision is given by 

m 5 b G 2 F \K tb K:fa 



T(B - X al ) 



32tt 4 



C 7 (fi b 



+ 



(8) 



where the auxiliary quantities CV(/x&) an d C7 (/-*&) are defined as 



C 7 (fi b ) = Cf(fx b ) 



C 7 b,g{^b) + —C7g,g(Hb) 
b 



c'M = cm - c 7Kg ^ b ) + —c 7 ,^b) 

ifib 



16V2n 3 a s ({ib) 
Gf K tb Kt s 



(9) 



where 



(10) 



Note that we have neglected the small contributions from the operators 7c ~ g and 0' 7c g . The 
branching ratio can be expressed as 



BR(5 -> X sl ) = r(E 7 Xs7) BR 5 l 



SL 



(11) 



where BRsx = (10.49±0.46)% is the measured semileptonic branching ratio. To the relevant 
order in a s , the semileptonic decay width is given by: 



ml Gl\K r 



SL 



L b^F \ xl cb\ / ™> c 



192tt 3 g \m 2 b/ 



(12) 



where the phase-space function g(z) is g(z) = 1 — 8z + 8r 3 — z A — Viz 2 log z. 

Note that there is no systematic distinction between the pole mass rrib and the corre- 
sponding running mass normalized at the scale \ib in the LL approximation. To be specific, 
the mass parameters are always treated as pole masses in our numerical analysis. 



III. SQUARK MASS MATRICES AS NEW SOURCES OF FLAVOUR 

VIOLATION 



As advocated in the introduction, the aim of this paper is to provide a phenomenological 
analysis of the constraints on the flavour violating parameters in supersymmetric models 
with the most general soft terms in the squark mass matrices. As explained there, we work 
in the mass eigenstate formalism, which remains valid (in contrast to the mass insertion 
approximation) when the intergenerational mixing elements are not small. 
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A specification of the squark mass matrices usually starts in the super-CKM basis, in 
which the superfields are rotated in such a way that the mass matrices of the quark fields 
are diagonal. In this basis, the 6x6 squared-mass matrix for the d-type squarks has the 
form 



f m 2 



Ml 



+ F dLL + D. 



d, LL f r dLL 



dLL 



{™\lr) +F dL R 



m. 



d, lr ) +FdRL rn dRR + F dRR + D dRR 



(13) 



For the it-type squarks we have 



M: 



m 



it, LL 



+ F uLL + D u LL (m 2 ^ LR )+F U LR 



, LR 



RL 



,RR 



+ F U 



RR 



D u 



(14) 



RR 



In this basis, the F terms (stemming from the superpotential) in the 6x6 mass matrices 
M. 2 (/ = u, d) are diagonal 3x3 submatrices, reading (Ff RL = Fj LR ) 



(F d LR)ij = ~K m d,i tan/3) % , {F uLR ) i:j = -//(m Uji cot (3) t 



(15) 



(F dLL )ij — m di tij , (F uLL )ij — m 2 ui \j . 



(16) 



Also the D-term contributions -D/ll and Df RR to the squared-mass matrix are diagonal in 
flavour space: 



D f ll,rr = cos 2(3 m| (t| - Q f sin 2 9 W ) 1 3 . 



(17) 



Since present collider limits give indications that the squark masses are larger than those 
of the corresponding quarks, the largest entries in the squark mass matrices squared must 
come from the soft potential, directly linked to the mechanism of supersymmetry breaking. 
These contributions, denoted in ( |T3"D and (0) by rn 2 j LLl rn 2 j RR and rn 2 f LRl are in general 
not diagonal in the super-CKM basis. 

Further comments are in order. Because of SU(2)i invariance, m 2 u LL and m 2 d LL are 
related. In the super-CKM basis this relation reads m 2 u LL = Km d LL K\ where K denotes 
the CKM matrix. The off-diagonal 3x3 block matrix m 2 j LR equals A d v d for down-type and 
A*v u for up- type squarks (the two vacuum expectation values are chosen to be real). They 
arise from the trilinear terms in the soft potential, namely A d ^H d DiDj and A U ^H U UiUj. 
We stress that we do not assume the proportionality of these trilinear terms to the Yukawa 
couplings, as is done in the mSUGRA model. Furthermore, differently from m 2 j LL and 
m 2 ^ RR , the off-diagonal 3x3 matrix m 2 j LR is not hermitian. 

Because all neutral gaugino couplings are flavour diagonal in the super-CKM basis and 
the mixing in the charged gaugino coupling to quarks and squarks is governed by the con- 
ventional CKM matrix, the flavour change through squark mass mixing is parametrized by 
the off-diagonal elements of the soft terms m 2 f LL , m 2 f RR , m 2 f LR in this basis. 



S 



The diagonalization of the two 6x6 squark mass matrices M 2 d and J\A 2 U yields the 
eigenvalues ra| and m\ k (k = 1,...,6). The corresponding mass eigenstates, and d k 

(k = 1,...,6) are related to the fields in the super-CKM basis, uij, URj and d^, d,Rj, 
(j = 1, ...,3) as 

u l ,r = T\j LtR u , d L)R = Y ] DL R d, (18) 
where the four matrices Tul,r and Tdl,r are 6x3 mixing matrices. 

IV. PHENOMENOLOGICAL ANALYSIS 

Our phenomenological analysis is based on a complete LL QCD calculation within the 
unconstrained MSSM; it is done in two parts: 

• In the first part, we try to derive bounds on the off-diagonal elements of the squark 
mass matrices by switching on only one of these elements at a time. We include, 
however, all new physics contributions (chargino, neutralino, charged Higgs bosons, 
gluino) in the analysis. We show that only those parameters get stringently bounded 
by B — > X s 7, which can generate contributions to the five- dimensional gluino-induced 
dipole operators 7g ~ g and 0' 7gg . 

• In the second part of our analysis we investigate whether the bounds obtained in the 
first part remain stable, if all off-diagonal elements, which induce the decay B — > X s 7, 
are varied simultaneously. We anticipate that the bounds on the individual off-diagonal 
elements get lost, because in this case various combinations of off-diagonal elements 
can contribute (with opposite sign) to the Wilson coefficients of the five-dimensional 
dipole operators. In the scenarios we discuss below, it is, however, possible to constrain 
certain simple combinations of off-diagonal elements of the down squark mass matrices, 
provided tan/5 and ji are not very large. 

A. General comments 

In order to analyse the implications of B — > X s j on the flavour violating soft parameters 
in the squark mass matrices, we choose some specific scenarios that are characterized by the 
values of the parameters 

li, M H -, tan/5, M susy , m g . (19) 

We regard this as reasonable, because we expect that these input parameters, which are 
unrelated to flavour physics, will be fixed from flavour conserving observables in the next 
generations of high energy experiments (provided low energy SUSY exists). Note that the 
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common SUSY scale, M susy , fixes in our scenarios the general soft squark mass scale rriq (see 
eqs. ( pTj|22D ) and the first diagonal element of the chargino mass matrix M 2 (see eq. (|A4|) ) . 



The parameters are chosen as follows: for M susy we choose the three values M susy = 
300, 500, 1000 Ge^, while for tan/5 we use the values tan/3 = 2,10,30,50. For the gluino 
mass, characterized by x = m?/M s 2 usy , we take x = 0.3, 0.5, 1, 2. Unless otherwise stated, the 
\i parameter and the mass of the charged Higgs boson M H - are fixed to be fi — 300 GeV 
and M H - = 300 GeV. 

While in the first part of our analysis we set, following ref. f|], all diagonal soft entries 
in m 2 dLL , m 2 dRR , and m 2 u RR equal to the common soft squark mass scale m~, we relax this 
condition in the second part of our analysis. 

We point out that the present bound on the mass of the lightest neutral Higgs boson 
requires a non-vanishing mixing (?tj.„ 1,^)33 among the stop-squarks. For our choices of the 
parameters, the MSSM bound coincides with the bound on the SM Higgs boson. 

We note that there are two contributions to the stop-squark mixing, namely the 'soft' 
contribution (m^ LR ) 33 and the F-term contribution (—fj,m t cot/?). In a general uncon- 
strained MSSM, the soft contribution does not scale with m t . However, following common 
notation, we parametrize the stop mixing term in hA\ (see eq. (|bj|)) as 

X t m t = {m 2 uLR ) m - iim t cot/3. (20) 

We fix the stop mixing parameter X t such that the mass of the lightest Higgs boson is at 
least 115 GeV to assure that the present Higgs bound is fulfilled in our analysis. We use the 
program FeynHiggsFast ]18[ which determines the Higgs boson mass approximately taking 



into account the complete one- and two-loop QCD corrections, the effects of the running 
top mass and of the Yukawa term for the light Higgs boson. The input parameters of the 
FeynHiggsFast program, are tan (3, the diagonal entries of the stop and the sbotton squark 
mass matrix, M~ l , M~ r , M~ , and M~ , the stop and the sbotton mixing parameters, X t and 
X , the top mass m t , the parameter fi and the charged Higgs boson mass. For tan/3 = 10 
we choose the following values for X t in dependence of the parameter M susy : (M susy ,X t ) = 
(300 GeV, 470 GeV); (500 GeV, 750 GeV); (1000 GeV, 1200 GeV). (For fx = 300GeV and 
X = we find Higgs boson masses of 115.2 GeV, 119.9 GeV and 121.1 GeV , respectively.) 
The dependence of the Higgs boson mass on the parameters 11 or Xb is rather small within 
our parameter scenarios. Also for the tan/3 = 30 and 50 scenarios we use the same values 
of X t for our convenience. In these cases the chosen X t values imply slightly higher Higgs 
boson masses. We also note that within our choice of parameters, the low tan (3 scenario, 
with tan (3 = 1, already gets excluded by the bound on the Higgs boson mass. 



B. First part of analysis 



In this part of the analysis only one off-diagonal entry in the soft part of the squark 
mass matrices is different from zero. We further assume (as in ref. ||) that all diagonal soft 
entries in m 2 d LL , m 2 dRR , and m 2 u RR are set to be equal to the common value m? = M 2 usy . 
Then we normalize the off-diagonal elements to m?, 
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FIG. 1. Dependence of BR(B — ► X s y) on 8d,LR,23- I n the upper frame, only SM and gluino 
contributions are considered. In the lower frame, the additional contributions from chargino, 
charged Higgs boson and neutralino are included. The horizontal lines denote the experimental 
limits. The different lines correspond to different values of x = rrv~/m~: 0.3 (short-dashed line), 0.5 
(long-dashed line), 1 (solid line), and 2 (dot-dashed line). The other parameters are fx = 300 GeV, 
tan/3 = 10, M H - = 300 GeV , M snsy = 500 GeV. 



[m 



>f,LL,ij 



f, LLjij 



m\ 



s 



] f,RR,ij - 



_ ( m /, RRjij 



mi 



(i ^ h 



(21) 



[m 



7 f,LR,ij 



f, LRJij 



[m 



2 \t 
/, LRJij 



mi 



(22) 



We recall that the matrix m u ^ LL cannot be specified independently; SU{2) L gauge invariance 
implies that m u ^ = Km^^K^ , where K is the CKM matrix. We also note that our 8- 
quantities only include the soft parts of the matrix elements of the squark mass matrices, 
while in ref. |§ also the F-term contributions are included in the definition of the (^-quantities. 
In figs, [l] and ^|, we show the dependence of the branching ratio of B — >• X s 7 on the flavour- 
violating parameters 6d,LR,i3 and Sd,RL,23, respectively. The upper frame in each figure is 
borrowed from i.e. we consider only SM and gluino contributions. As 5a,LR,23 and 
8d,RL,23 generate the five- dimensional dipole operators O^-g-g and 0' 7gg , it is not surprising 
that they get stringently bounded. We should note that at this level of the analysis there 
is no dependence of these bounds on fi or tan (3. Such a dependence could result from the 
term (F dtLR ) 33 , but only when 5d,LL,23 and 5d,i?_R,23 are turned on. We will discuss this point 
in more detail at the end of this section and in the second part of our analysis. In the 
lower frame of figs. [I] and we also include the contributions from charginos, charged Higgs 
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FIG. 2. Same as in fig. [I] when 5d t RL,23 is the only non-zero off-diagonal squark mass entry. 



bosons and neutralinos. Comparing the branching ratio in the two frames at 5d,LR,23 — 
and 5 d,RL, 23 — (which corresponds to switching off the gluino contribution), one concludes 
that the combined contribution from charginos, neutralinos and charged Higgs bosons is of 
the same order as the SM contribution. A detailed investigation shows that the neutralino 
contribution is negligible, while the contributions from the chargino and charged Higgs 
boson are similar in magnitude; both interfere constructively with the SM contributions for 
the specific choice of parameters. However, as the gluino yields, intrinsically, the dominant 
contribution by far, the bounds 5d,LR,23 and Sd,RL,23 are only marginally modified by chargino, 
neutralino and charged Higgs boson contributions. A comment concerning the different 
shapes of the curves in figs. [I] and |2] is in order. In fig. |^, with non-vanishing 5d,RL,23i 
the gluino contribution is induced by the primed-type operator O'^- - and therefore does not 
interfere with the contributions from the other particles, as these induce unprimed operators 
in the first place. In contrast, in fig. |, which shows the case of non-zero 5d,LR,23, the gluino 
contribution is of the unprimed type and therefore interferes with the other contributions. 

We also tried to derive analogous bounds on 5 d ,LL,23, 8d,RR,23, $u,lr,23, 8 UtRL ,23, $u,rr,23 
and also on S U) lr^3 and S U) lr j2 2- In the chargino sector the latter diagonal elements, together 
with the usual CKM mechanism, also can induce flavour violation. The parameters of the up- 
squark mass matrix give rise to chargino contributions that lead only to dimension six dipole 
operators, which inherently are not very large. For our choices of /i, M susy and tan (5, this was 
confirmed numerically. Therefore, no stringent bounds are obtained for the soft parameters 
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in the up-squark mass matrix^. The remaining parameters of the down-squark mass matrix, 
i-e. Sd,LL,23 and S^rr^, play an interesting role. They not only generate contributions to the 
six-dimensional operators in (|6|), but, together with the chirality changing term (F djLR ) 33 , 
they also induce contributions to the five-dimensional gluino operators in (0). For the values 
of \x and tan j3 used in our analysis, the coefficients of the five-dimensional operators turn 
out to be rather small. Thus, no stringent bounds on 6d,LL,23 and 5d,RR,23 are obtained. 

Summarizing the first part of our analysis, we conclude that 5d,LR,23 and 5d,RL,23 are 
the only parameters that get significantly constrained by the measurement of the branching 
ratio of B — > X s 7. 



C. Second part of analysis 



We now explore the problem of whether the separate bounds on Sd t LR,23 and 5d,RL,23, ob- 
tained in the first part, remain stable if the various soft parameters are varied simultaneously. 
The analysis is based on the assumption that the soft terms in the squark mass matrices 
have the hierarchical structure that the diagonal entries in m 2 d LL , m 2 d RR , and m 2 u RR are 
larger than the off-diagonal matrix elements (including m 2 d LR and m 2 u LR ). In contrast to 
the first part of the analysis, we will allow for a non- degeneracy of the diagonal elements 
in the matrices m 2 dLL , m 2 dRR , and m 2 u RR . To implement this, we define (^-quantities in 
addition to those in eqs. (f2Tj) and (|2"2"|), which parametrize this non-degeneracy: 

x {m 2 LLL )u-m 2 g {m\ RR )n-m 2 g 

d f,LL,ii ~ ^2 ' °f,RR,ii ~ ^2 ' \ Z6 > 

Unless otherwise stated, the diagonal (^-parameters (in eq. (|23|)) are varied in the interval 
[—0.2, 0.2]. On the other hand, the off-diagonal ones (in eqs. (^) and (p2|)) are varied in the 
interval [—0.5,0.5], by use of a Monte Carlo program. There are, however, two exceptions. 
First, we do not vary those off-diagonal <5's with an index 1; the latter <5's we set to zero, since 
they are severely constrained by kaon decays (see for example ||). Second, as mentioned 
earlier, also (m^ LR )33 is not varied, but fixed such that the mass of the lightest neutral Higgs 
boson gets heavy enough to be compatible with experimental bounds. 

In our Monte Carlo analysis we plot those events, corresponding to 2.0 x 10~ 4 < BR(B — > 
X s >y) < 4.5 x 10" 4 , which is the range allowed by the CLEO measurement. Note that we 
do not include recent preliminary data || in our analysis. Furthermore, we have made sure 
that our events correspond to squark masses that are real and lie above 150 GeV. The 
dependence of the bounds on this specific choice is discussed below. We start with the 
following parameter set: \i = 300 GeV, M H - = 300 GeV, tan/? = 10, M susy = 500 GeV, 
x = m? / M 2 usy = 1 and X t = 750 GeV. In fig. [3|, we only consider SM and gluino 



x In [19] the authors derived a rather stringent bound on a quantity proportional to <5 Ui Li?,,33 in 
the case of a small chargino mass of 100 GeV. However, they include the small CKM factor 



K^Kfb ~ 1/30 in the definition of their quantity. 
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FIG. 3. Contours in the 5d,LR,23 - <>d,RL,23 plane. In the left frame, 5d,LR,23 and 5d,RL,23 are the 
only flavour-violating parameters. In the right frame, we allow also for non-vanishing Sd n 23 and 
Sd,RR,23- We only consider SM and gluino contributions and the other parameters are fj, = 300 GeV, 
tan/3 = 10, M susy = 500 GeV and x = 1. 



contributions. In the left frame we present the constraints on 6d,LR,23 and 5d,RL,23 when 
these are the only flavour-violating soft parameters; the diagonal (^-parameters defined in 
eq. fl23|) are also switched off. As expected from the first part of our analysis, stringent 
constraints are obtained. The hole inside the dotted area represents values of 5d,LR,23 and 
5d,RL,23 for which the branching ratio is too small to be compatible with the measurements. 
In the right frame we investigate interference effects that arise when 6d,LL,23 an d $d,RR,23 are 
switched on in addition to 5d,LR,23i &d,RL,23- All of them are varied between ±0.5. From fig. |3| 
we find that the bounds on 5d,LR,23 and 5d,RL,23 cannot be softened significantly by non-zero 
values of 5d,LL,23 and Sd,RR,23, although these <5-parameters, which individually give rise to 
six- dimensional operators, generate five- dimensional operators through the interplay with 
the F-term (F d ,LR) 33- As already discussed in the first part of the analysis, for moderate 
values of /1 and tan/3, the contribution to the Wilson coefficient of the five- dimensional 
operator is rather small. 

The full power of the interference effects from different sources of flavour violation is 
depicted in fig. g where we allow not only for non-zero 5 d ,LR,23, &d,RL,23, S d ,LL,23 and 5d,RR,23 
but also for non- vanishing 5d,LR,22, 8d,LR,33- All these parameters are varied between ±0.5. 
As can be seen, the bounds on 5d,LR,23 and 8d,RL,23 get destroyed dramatically. The reason 
is that there are now new contributions to the five- dimensional dipole operators. As an 
example, the combined effect of 6d,LR,33 and 6d,LL,23 leads to a contribution to the Wilson 
coefficient of the operator 0-jg~ g . The sign of this contribution can be different from the one 
generated by &d,LR,23- As a consequence, the bound on 5d,LR,23 gets weakened. To illustrate 
this more quantitatively, we assume for the moment that there are only these two sources 
that can generate 07§,§, i.e. we switch off the other 5- quantities. If 5d,LR,23 is larger than 
the individual bound from the first part of the analysis, it is necessary that the product of 
^d,L_R,33 and Sa,LL,23 1S also relatively large; only in this case can the two sources lead to a 
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FIG. 4. Contours in the S djLR) 23 - S d)RL , 2 3 plane, where 5d,LR,23, S d ,RL,23, &d,LL,23 and 5 dtRR}2 3 
are the only flavour- violating parameters; 5d,LR,22-, $d,LR,33 are also non- vanishing. We only include 
SM and gluino contributions; the other parameters are \i = 300 GeV , tan = 10, Af susy = 500 GeV 
and x = 1. 



branching ratio compatible with experiment. This feature is illustrated in fig. only values 
of S d 

l R ,23 an d values of Sd,LR,33 ' $d,LL,23 which are strongly correlated lead to an acceptable 
branching ratio. As clearly visible from fig. |5|, the correlation between the two sources for 
Org,g, is essentially linear. This implies that the linear combination 

$d,LR,23 + f$d,LR,33 " $d,LL,23 (24) 

gets constrained, if / is chosen appropriately. Stated differently, the Wilson coefficient of the 
operator Or§,g is essentially proportional to the combination fl24f) . This implies in turn that 
for the values of the parameters we are using at the moment (/i = 300 GeV , M H - = 300 GeV, 
tan/5 = 10 M susy = 500 GeV, x = m? / M s 2 usy = 1, X t = 750 GeV), the Wilson coefficient is 
well approximated by its double mass insertion expression. The coefficient /, which can be 
read off from this expression, depends on the parameter x = m?/M s 2 usy and reads 

1 + 9x — 9x 2 — x 3 + 6x(l + x) logx 
~ (1 -x)[5x 2 -Ax- 1 -2x(x + 2)logx] ' ^ ' 

The numerical values of f(x) for some values of x read 0.74 for x = 0.3, 0.68 for x = 0.5, 
0.60 for x = 1.0 and 0.52 for x = 2.0, respectively. 

The solid line in fig. || represents pairs (Sd,LR,23, ^d,LR,33 • $d,LL, 23) for which the combi- 
nation in eq. (£23|) is zero. The points scattered around this line therefore represent Monte 
Carlo events for which this combination is small. We now turn back to the scenario of fig. 
I in which all the parameters 5 d ,LR,23, S d>RL ,23, <$d,LL,23, ^,^,23, $d,LR,22, S d ,LR,33 are varied 
simultaneously. In this case, the linear combinations 

LCi = 8d,RL,23 + f( x )3d,RR,23 " $d,RL,33 + f ( x )$d,RL,22 ' 8d,LL,23, 

LC 2 = S d ,LR,23 + f(x)5 d ,LR,22 " S dlRR ,23 + f(x)S d ,LL,23 " $d,LR,33, (26) 



15 




d,LR,23 



FIG. 5. The parameters 5d t LR,23 and 5d,LR,33 • $d,LL, 23, which are compatible with the data 
on B — > X s 7 are shown by dots. Values lying on the solid line lead to a vanishing contribution 
of the five-dimensional operator O^gg in the MIA. See text. We only switch on SM and gluino 
contributions; the other parameters are fi = 300 GeV, tan/3 = 10, M susy = 500 GeV and x = 1. 



are expected to get constrained. 

In fig. ^|we show the allowed region for LG\ and LG 2 . There, we allow all non-diagonal 
5-parameters to vary between ±0.5. In addition, we also allow for non-equal diagonal soft 
entries, by varying the parameters 5f,LL,u and S^rr^h between ±0.2. With the latter choice 
we still guarantee the hierarchy between diagonal and off-diagonal entries, but we get rid of 
the unnatural assumption of degenerate diagonal entries. In the left frame, we include only 
SM and gluino contributions. We find that the linear combinations LC\ and LC 2 indeed 
get stringently bounded. In the right frame of fig. ^| we test the resistance of these bounds 
when the additional contributions (i.e., those from charginos, charged Higgs bosons and 
neutralinos) are turned on. In this case also 5 UiLRi2 3, S U)RLj23 , 5 UiLLi23 , S UjRR>23 and S u , L r,22 
are varied in the range ±0.5. We find that the bound on LC\ remains unchanged, while 
the one on LC 2 gets somewhat weakened. This feature is expected, because charginos and 
charged Higgs bosons contribute to unprimed operators at first place. At this point we should 
stress that these plots were obtained by choosing the renormalization scale fib = 4.8 GeV 
and by requiring all squark masses to be larger than 150 GeV. We checked that the bounds 
on LC\ and LC2 remain practically unchanged when the renormalization scale is varied 
between 2.4 GeV and 9.6 GeV; they are also insensitive to the value of the required minimal 
squark mass, as we found by changing m squar k m in from 150 GeV to 100 GeV or 250 GeV. 
Moreover, we also checked whether the restriction to fi = ±300 GeV scenario is too severe: 
we redid the complete analysis for \i = —300 GeV and confirmed that ther are no differences 
between the results of these two choices. 

Two remarks are in order: First, one might wonder why we did not include terms like 
$d,RR,33 ■ 8d,LR,23 m LC\ and LC 2 , which would result into more complicated combinations. 
As we are allowing for nonequal diagonal soft entries, these terms give in principle additional 
contributions to the five dimensional operators. However, as the diagonal ^-parameters are 
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FIG. 6. Contours in the LCi-LC 2 plane with 5 d;LRi23 , S d;RLi23 , 5 d ,LL,23, S diRR ,23, S diLR)22 , 
$d,LR,33, 8 u ,lr,23, &u,RL,23, S u ,ll,23, S U)RR)23 , and 5 U)LRi22 all non- vanishing. In the left frame, 
we consider only SM and gluino contributions whereas in the right frame we also include chargino, 
charged Higgs boson and neutralino contributions. The values of the other parameters are 
fi = 300 GeF, tan/3 = 10, M H - = 300 GeV , M susy = 500 GeV and x = 1. 



only varied between ±0.2, their influence on the Wilson coefficients is numerically small. 
For this reason, the simpler combinations LC\ and LC 2 , denned in eqs. fl26|) , are sufficiently 
constrained and we prefer to give bounds on these quantities. 

Second, if we got rid of the hierarchy of diagonal and off-diagonal entries in the squark mass 
matrices, stringent bounds on the simple combinations LC\ and LC 2 certainly would no 
longer exist, simply because there would then be more contributions to the five- dimensional 
operators of similar magnitude. In this case, however, the full Wilson coefficients of the 
five-dimensional operators still would be stringently constrained by the experimental data 
on B — > Xg'y. Unfortunately, in this case not much information can be extracted for the 
individual soft parameters or simple combinations thereof. 

Finally, we extend our analysis to other values of the input parameters. So far, we 
found that the combinations LC\ and LC 2 (see eqs. (|26|)) are stringently bounded in the 
scenario characterized by the input values \i = 300 GeV, M H - = 300 GeV, tan/3 = 10, 
M susy = 500 GeV, x = m| / M s 2 usy = 1 and X t = 750 GeV. It is conceivable that the 
bounds on LC\ and LC 2 can get considerably weakened in other scenarios. Therefore, we 
analyse the bounds on the soft parameters within the following parameter sets: (M susy , X t ) 
= (300 GeV, 470 GeV), (500 GeV, 750 GeV), (1000 GeV, 1200 GeV). For tan/3 we explore 
the values: tan/3 = 10, 30, 50. Furthermore, the gluino mass nig is varied over the values 
x = m?/M s 2 usy = 0.3,0.5,l,2. 

Surprisingly, the constraints on LC\ and LC 2 are completely stable over large parts of 
the parameter space. Within the tan /3 — 10 scenario the bounds are essentially unchanged 
if the other two parameters M susy and x, are varied over the complete range of values given 
above. For example, the independence from the parameter M susy within this scenario can 
be read off from the comparison of frames in the first vertical line in fig. |7]. 
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FIG. 7. LCi versus LC2 in the x = 1 scenario. Vertically we vary the parameter M susy from 
M susy = 1000 GeV in the first line, via M susy = 500 GeV in the second to M susy = 300 GeV in 
the last. Horizontally we vary tan/3 from 10 via 30 to 50. The values of the other parameters are 
/i = 300 GeV and M H - = 300 GeV . All contributions are switched on. 



However, fig. |7] also illustrates that the bounds get significantly weakened or even lost 
when tan /3 values as large as 30 (second vertical line) or 50 (third vertical line) are chosen. 
This effect gets enhanced when the general mass scale m§ in the squark mass matrices 
decreases with the parameter M susy . There are two main reasons why the bounds get 
weakened in these scenarios. First, in the large tan/3 regime the term (Fd,LR)33 gets strongly 
enhanced because of its proportionality to tan/3 (see ([ToD). Particularly, for tan/3 = 50 
and M susy = 300 GeV, the term (F dt L R ) 33 is of the same magnitude as the diagonal entries 
of the squark mass matrix. Thus, the contributions to the Wilson coefficients of the five- 
dimensional gluino operators (induced by (F^lr) 33 in combination with 5d,LL,23 or 8d,RR,2s) 
become important enough to weaken the bounds on LG\ and LC2 significantly. The relative 
importance of this F term is of course increased if the general soft squark mass scale M susy 
is decreased as can be read off from fig. [7|. Second, within the large tan/3 regime the 
contributions from charginos get enhanced and therefore also weaken the bounds on LCi- 
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FIG. 8. The left frame shows the bounds on LC\ and LC2 in a tan (3 = 50 scenario, which are 
relatively weak. In the right frame we put (F c i ! lr)33 = 0. See text. The other parameters are 
/i = 300 GeF, M H - = 300 GeV, M snsy = 500 GeV and x = 1. All contributions are switched on. 



These features are illustrated in more detail in fig. || In the first frame we take over the 
specific scenario with tan j3 = 50 and M susy = 500 GeV from fig. [7|. To show that the term 
{Fd,LR) 33 is indeed one of the reasons for the weakening of the bounds, we present in the 
right frame of fig. ||the corresponding scenario when (Fd t LR)33 is set to zero. We see that we 
regain better bounds on LC\ and also on LC2. However, we also see that the bound on LC2 
remains weak. This, and the resulting asymmetry, is due to a large chargino contribution for 
tan j3 = 50. We recall that there is no chargino contribution to the primed operator which 
could influence the bound on LC\. 

We can also explore how the bounds behave if we vary the parameter [i. Until now we 
used the value fi = 300 GeV. Because the parameter (F djLR ) 33 is actually proportional to 
the product of tan/3 and /i (see eq. ([l^)), we conclude from the findings above that the 
bound on LC\ is unchanged if we increase the value of \i and decrease the value of tan (3 such 
that the product of both parameters is constant; the bound on LC2 is then even stronger 
because the chargino contribution is smaller for increasing fi. Consequently, one finds a 
smaller asymmetry in the corresponding plots (compare the left frame in fig. ^ with the 
second frame in the second line of fig. 0). On the contrary, if one decreases the value of \i 
to fi = 150 GeV, the bound on LC2 is weakened and the asymmetry of the plot is increased 
as one can read off from the right frame in fig. 

Summing up the second part of our analysis, the two simple combinations LC\ and 
LC2 (Prj|), consisting of elements of the soft parts of the down squark mass matrices, stay 
stringently bounded over large parts of the supersymmetric parameter space, excluding the 
large tan (3 and the large [i regime. We note that these new bounds are in general one order 
of magnitude weaker than the bound on the single off-diagonal element Sd t LR,23, which was 
derived in previous work P,pO| by neglecting any kind of interference effects (see e.g. tab. 4 



in [20J where the value 1.6 TO 2 is given as bound on ^^^23 for x = 1 and M susy = 500 GeV). 
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FIG. 9. The left frame shows bounds on LC\ and LCi for a scenario with /i = 900 GeV 
and tan/3 = 10. The right frame shows the corresponding bounds obtained with fi = 150 GeV 
and tan/3 = 30. The other parameters in both frames are chosen to be M H - = 300 GeV, 
-^susv = 500 GeV and x = 1. All contributions are switched on. See text. 



V. IMPLICATIONS ON C 8 (fi w ) AND C 8 (/i w ) 



As mentioned in section |J it is possible to absorb the various versions of gluonic dipole 
operators into the SM operator 0$ and its primed counterpart. The resulting effective 
Wilson coefficients, denoted by C$(fj,w) and C' 8 (nw), read at the matching scale fiw'- 

G8{fJ-W) - ^8{^W> - ^8b,g{^W) H j T C ' 8 g ,g{^W ) ^ 77 ^ 

\ m b {/j, w ) J LxF&tb&ts 

C's(nw) = C&iw) ~ (cfaM + Cwfrw)) l ^ K as ^ w) . (27) 

\ ' y m b {iJL W ) J G F K th K* ts 

The coefficients on the r.h.s. of eq. ([27D are given explicitly in section [A 2| (appendix 

We now investigate the implications on possible values for the effective Wilson coefficients 
C%(nw) and C' 8 (/j,w) when taking into account the experimental constraints on B — > X s ^y. 
The result is shown in fig. 0, for \i = 300 GeV, M H - = 300 GeV, M susy = 500 GeV, 
X t = 750 GeV, tan/3 = 10 and x = 1. The soft parameters, encoded in the 5 quantities, 
are varied as in fig. |6|. From fig. |10| we conclude that large deviations from the SM values 
for Cs(fiw) and C' 8 (fiw) are still possible. Scenarios in which these Wilson coefficients are 
enhanced with respect to the SM gained a lot of attention in the last years. For a long time 
the theoretical predictions for both, the inclusive semileptonic branching ratio B s i and the 
charm multiplicity n c in i?-meson decays were considerably higher than the experimental 
values [T2T[] . An attractive hypothesis, which would move the theoretical predictions for both 
observables into the direction favoured by the experiments, assumed the Wilson coefficients 
Cs(nw) and C' 8 (/j,w) to be enhanced by new physics W2 . 
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FIG. 10. Allowed values for the Wilson coefficients C&(my/) and C' s (mw) satisfying the ex- 
perimental constraint on BR(B —* X s j). The values of the input parameters are \i = 300 GeV, 
tan/3 = 10, M H - = 300 GeV, M susy = 500 GeV and x = 1. All contributions are switched on. 



After the inclusion of the complete NLL corrections to the decay modes b —>■ cuq and 
b — > ccq (q = d, s) p3| , the theoretical prediction for the central values of the semileptonic 
branching ratio and the charm multiplicity are still somewhat higher than the present 
, but theory and experiment are in agreement within the errors. It should 
however, that in the theoretical error estimate the renormalization was varied 



measurements 
be stressed, 

down to rrifj/A. If one only considers the variations down to the theoretical predictions 

will have only an marginal overlap with the data. This implies that there is still room for 
enhanced C 8 ([J>w) an d C' 8 (fiw) PU| . 



VI. SUMMARY 



In this paper we have chosen the rare decay B — > X s 7 to analyse the importance of 
interference effects for the bounds on the parameters in the squark mass matrices within the 
unconstrained MSSM. Our analysis, based on a systematic leading logarithmic (LL) QCD 
analysis, mainly explored the interplay between the various sources of flavour violation 
and the interference effects of SM, gluino, chargino, neutralino, and charged Higgs boson 
contributions. Surprisingly, such an analysis did not exist so far. Unlike previous work, 
which used the mass insertion approximation, we used in our analysis the mass eigenstate 
formalism, which remains valid even when some of the intergenerational mixing elements 
are large. 

In former analyses no correlations between the different sources of flavour violation were 
taken into account. Following that approach, we found only two down-type squark mass 
entries to be significantly constrained by the data onB-> X s ^\ 5d,LR,23 and 5d,RL,23- These 
entries are correlated with the five- dimensional dipole operators where the chirality flip is 
induced by the gluino mass. We showed that these bounds get destroyed in scenarios in which 
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certain off-diagonal elements of the squark mass matrices are switched on simultaneously. 

We then systematically explored the interference effects from all possible contributions 
and sources of flavour violation within the unconstrained MSSM. Accordingly, we switched 
on all off-diagonal elements 5i of the squark mass matrices and varied them in the range 
±0.5. In addition, we also varied the diagonal elements, but in smaller interval in order to 
preserve a certain hierarchy between the off-diagonal and the diagonal ones. In this general 
scenario we singled out two simple combinations of elements of the soft part of the down 
squark mass matrix, which stay stringently bounded over large parts of the supersymmetric 
parameter space, excluding the large tan f3 and the large \i regime. These new bounds are 
in general one order of magnitude weaker than the bound on the single off-diagonal element 
^d,L_R,23, which was derived in previous work 0,^] by neglecting any kind of interference 
effects. 

Finally, we briefly analysed up to which values SUSY contributions, compatible with 
B — > X s 7, can enhance the Wilson coefficients Cs(mw) and C' s (mw)- We found that large 
deviations from the SM values are still possible in our general setting. Such scenarios are of 
particular interest within the phenomenology of inclusive charmless hadronic B decays. 
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APPENDIX A: 



1. Mixing matrices, interacting Lagrangian 



In this appendix we present our conventions in the mass mixing matrices for the relevant 
particles and in the interacting Lagrangian. Besides some specific changes, we follow |0j] : 

Charged Higgs bosons: If we denote the two SU (2) Higgs boson doublets appearing in 
the superpotential by 

^=(3)- < ai > 



the corresponding mass eigenstates 
by (see |fT6| ) 



and B.2 of the charged Higgs bosons are given 



Hi* 
Hi 



sin (3 cos (3 
- cos (3 sin (3 




(A2) 



Ze 
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and similarly for H 2 = (H 2 )* an d = (Hi)*. 

In the unitary (physical) gauge, the massless charged fields H 2 ~ are absorbed by the W 
boson. One is left with two massive charged Higgs bosons of equal mass. 

Charginos: The charginos xt/2 are a mixture of charged gauginos and Higgsinos K{ 
and h 2 . Defining 





(A3) 



the mass terms are then Ct = —^(i/j +t X t i/j + ip T Xil) + ) + h.c, where 

x = ( M 2 g 2 v 2 \ 
V 92V1 n J ' 



(A4) 



The two-component charginos xt (i — 1? 2) and the four- component charginos xt/2 are then 
defined as 



Xi 



xt 



(- 

V xT 



, xt 




(A5) 



where the unitary matrices U and V diagonalize X: Mf} — U*XV 1 = VX^U T . then 

becomes £^ = — M$ n xtxt ~ -^022 xtxt • U an d ^ can be found by observing that 
M% 2 = VX T XV- 1 = U*XX T U*- 1 . They are not fixed completely by these conditions. 
The freedom can be used to arrange the elements of Mf} to be positive: If the i th eigenvalue 
of is negative, simply multiply the i th row of V with —1. 

Neutralinos: The neutralinos are linear combinations of the gauginos A' and A 3 and the 
neutral Higgsinos h\ and h 2 . If we define 



( -iA' \ 
-i\ 3 

K 

\ hi J 



(A6) 



the neutralino mass term reads = —\ip 0T Yip Q + h.c, where 



Y 



I 


Mi 





givi 


91V2 
V2 


\ 







M 2 


g2v\ 

V2 


92 V 2 
V2 






gx vi 


92 vi 
V2 





-fi 




\ 


giV2 

V2 


92 V 2 
V2 


-n 





/ 



(A7) 



Two- and four-component neutralinos must be defined as 
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X? = JV«V?, (i = l,...,4) 



A 7 



(A8) 



To diagonalize the mass matrix, A" must obey Njj = N*YN~ X , where Nd is a diagonal ma- 
trix. A" can be found, using the property N D = NY^Y N~ x . The eigenvalues and eigenvectors 
are found numerically. Possible negative entries in Nd are turned positive by multiplying 
the corresponding row of A^ by a factor of i. 

Quarks: The situation in the quark sector is in almost complete analogy to that of the 
SM. The quarks get their masses from the Yukawa potential when the Higgs bosons acquire 
a vacuum expectation value. We define the mass eigenstates by 



(A9) 



(AlO) 



(All) 



^li = U t L jU Lj , u { m = UgiiRj, 

a Li — u ij a Lj-, «ffi — u ij a Ry 

The mixing matrices must satisfy (i = 1,2,3) 

D R\dT D L] = \ D = diagf^f 

jjRyuTjjL] = \ D = diag(^ 

where 

v 2 = v^^sin/3. 

As can be seen, the eigenvalues of A" and X d are fixed by the quark masses and the minimum 
of the Higgs potential. In the SM, the only observable effect of the mixing is encoded in the 
CKM matrix K = U L D L \ appearing in the charged current. Therefore it is possible and 
convenient to set D L = D R = U R = 1L (=>■ U L = K). To be more precise, \ d and A e are 
chosen to be diagonal and A" = diag {^f^j K T . Although in our theory the mixing matrices 
appear in all kinds of combinations, we adopt this convention here, emphasizing that it is a 
choice made just for convenience. An underlying theory should fix the values of A" and \ d 
at some (high) scale. Note that in the main text we neglect the superscript m for the mass 
eigenstates. 

Squarks: If supersymmetry were not broken, squarks would be rotated to their mass basis 
with the help of the same matrices as their fermionic partners. In a more realistic setting 
we need to introduce a further set of unitary rotation matrices. The notation must be set 
up carefully because the mass eigenstates of squarks and sleptons are linear combinations 
of the partners of left- and right-handed partners of the corresponding fermions. The exact 
form of the mass matrices and the notation for the corresponding diagonalization matrices 



can be found in section III. 



Interaction Lagrangian: In order to fix further conventions we quote the relevant parts 
of the interaction Lagrangian: 
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Charged Higgs boson-quark-quark 

£ qq H — {^i>)rne.K mk diP L Uk{Z E )2nH i + ( \%) m eK^ k Uk~P R de(Z E )l n {H i )* 

+ (\ u D ) me u- e P L d m (Z E y in (Hry + \ u D ) me d~P R u e (Z E ) ln Hr. (A12) 



Note that in our basis, the terms proportional to the \ d D always come together with 
the CKM matrix K, while the A^ terms do not. 



S quark- quark- chargino 



qqx ch 



Ujdi 



A-ijePi + B? je P R 



Xt c + u\xt c A^Pr + B^Pl dj 



(A13) 



( \ d Ft "\ TT* Ud 



where Af je - V /v D ± V LJij w £2i ^iji 
and xf 1 ° denotes the charge-conjugated field. 



Squark-quark-neutralino 



£qqX° 



-djdi 



d~iXe Cf},P R + D jAf Pj 



ijt r L 



d 



J ■ 



1(1 T7 5 ) 



(A14) 



where = {\ d D T ] DL ) %J N* n - V^^r^iY^ 
A* = (^DRhN i3 +^ DLij ((2Q d + l) 9l N ei - g 2 N e2 ). 



2. Wilson coefficients 



We recall the Wilson coefficients at the matching scale /j,w- The non- vanishing Wilson 
coefficients for the SM are, at leading order in a s (x tw = mf/m^r) : 

C 2 sm{^w) = 1 

Cjsm^w) = 24 ^ tW _ x y (~ 8x tw + 3a; L + 12x tw - 7 + (lSx 2 w - I2x tw ) \nx tw ) 

C 8 sm{^w) = g (J tW _ x~)4 (~ x tw + ® x tw ~ 3x tw - 2 - 6x tw Inxtw) ■ (A15) 

The contributions from charginos, neutralinos and charged Higgs bosons match onto the 
(chromo)magnetic operators of the SM and the corresponding primed operators, which differ 
from the SM ones only by their chirality structure. The corresponding Wilson coefficients 
become somewhat involved |T5[ as they include many mixing matrices, whose definitions 
were given in appendix Al. One gets (using the abbreviation V = {4G F K tb KQ / 

C-j{nw) = Cjsm{^w) 

- - [cot 2 (3 xtH{QuFi{x tH ) +F 2 (x tH )) + 
XtH(QuF 3 (x tH ) + F±(x tH ))] + 
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+ 
+ 
+ 



1 1 



d r>d* 



T3 a R 



2Vm\ 2je 3je 
1 l m ch 

Xi r)d Ad* 

2Vm\. m b *W A W 



F i( x x f Uj ) + Q u F 2 (x 



X(„. 



F 3 (x x ch u .) + Q u Fi(x y ch u . 



XT u i> 



QaJ_ 

2Vm 2 ~ 

dj 



D d 2]i D^F 2 {x x% ) + ^D^C^F^) 



(A16) 



+ 



+ 



- [cot 2 (3 x m Fi(xtH) + x tH F 3 (x t H)] 
1 1 



2V m\ 
1 1 



B 2 jtB^ t F 2 (x x ch u .) + — - B 2je A 3 * e F 4: (x x ch l 



D%iDt 3 iF 2 {x x » dj ) + ^D d 2je C d 3 * e F 4 (x x o dj ) 



C 7 (/j,w) — — t: 



"' ' ' ' " b tan 2 (3 x tH {Q U F 1 (x tH ) + F 2 (x tH ) ) 



+ 



2 to 2 
1 1 



2Vmi. 



—A 2jl Af jt Fi(x x ch Uj ) + Q u F 2 (x x ch u .) 



X\ Uj , 



1 l m ch 

, X-t Ad r)d* 



+ 



2Vml m b 
Qd 1 



2Vm 2 r 



yd nrf* 



C 2ji C 3 * e F 2 (x x o d .) + —^ C 2ji D 3 )iF 4 {x x ^ d . 



1 m s m h ^ 2 

2 TO? 



tan 2 x tH Fi{x tH ) 



+ 



+ 



1 1 

2l7mi 

1 1 

2l7mT 



T7W ch 

A 2jZ A 3*jt F 2 {x x ch Uj ) + —^-A 2je B 3 * e F 4 (x x ch u .) 



C d 2ji C% e F 2 (x x o dj ) + ^C d 2jt Df jE F,{x x% ) 

1 1 L b 



(A17) 



(A18) 



(A19) 



where Q u = 2/3 and Q d = —1/3. We kept the charged Higgs boson contribution to the 
primed operators since they are proportional to tan 2 /? which could compensate the m s /m t 
suppression. The functions F^x) are defined at the end of this section. Although the Wilson 
coefficients C' 7 (fib) and C'^fib) of the primed operators are usually small, we retain them in 
our analysis. 

Among the coefficients arising from the virtual exchange of a gluino, the most important 
ones are those associated with the (chromo)magnetic operators: 

C 7b ,g(^w) = -Jq^ g ( T DL r DL) F 2 {x gdk ) , 

k=l d k 

C 7 g,g(^w) = m~ g - J2 ( t dr t *dl) F 4 (x g d k ) , 

k=l d k 
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C%b,g{l^w) 



167T 2 
1 



1 

E__^ ( -pkb -p* ks\ 
~~2 \ DL DL ) 



m- 



■^F 2 (x gdh ) - -F x (x gdk ) 



C 8 g,g(fJ>w) ~ m g 



fc=l o! fc 



2 V 1 DR l DL 



1 

6 



(A20) 



Note that the coefficients Cjg^nw) an d Css.sC/^vk) are of higher dimensionality to compen- 
sate the lower dimensionality of the corresponding operators. The ratios x g d k are defined as 
x g d k = m|/m| . The Wilson coefficients of the corresponding primed operators (which are 

not small numerically) are obtained through the interchange T l ^ R <-> in eqs. ( |A20| ) . For 
the Wilson coefficients of the scalar /tensorial four-quark operators we refer to [ITJ]. Finally, 
we define the functions Fi appearing in the Wilson coefficients listed above: 



F x (x) 
F 2 (x) 
F 3 (x) 
FAx) 



12 (x 



12 (x 



— — (x 3 — 6x 2 + 3x + 2 + 6x log xj 
] -— (2x 3 + 3x 2 - 6x + 1 - 6a; 2 log . 



2 (x - 1) 



x 2 — Ax + 3 + 2 log x ) , 



2(x 



— — (x 2 — 1 — 2x log x 



(A21) 
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